Quantum dot as thermal rectifier 
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We report the observation of thermal rectification in a semiconductor quantum dot, as inferred 
from the asymmetric hne shape of the thermopower oscihations. The asymmetry is observed at 
high in-plane magnetic fields and caused by the presence of a high orbital momentum state in the dot. 
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Recently, the field of thermoelectricity and solid state 
thermionics has gained renewed attention, due to ad- 
vances in growth and fabrication of complex compounds, 
mesoscopic devices, and nanostructures [1, 2]. The main 
idea is to enhance the efficiency of macroscopic devices 
by the control of the energy transport on a microscopic 
scale. For example, a recent proposal to build a thermal 
rectifier in a non-linear lattice attracted much attention 
[3] . An actual nanoscale solid state thermal rectifier using 
tailored carbon and boron nitride nanotubes was recently 
demonstrated [4]. Here, we present yet another way to 
obtain thermal rectification, utilizing the thermoelectric 
properties of a semiconductor quantum dot (QD). 

QDs are the smallest possible thermoelectric devices. 
In the Coulomb blockade (CB) regime, their transport 
properties are highly nonlinear and strong thermoelectric 
signals result [5]. The thermopower 5'qd parametrizes 
the electric response Vt of a QD to an applied tem- 
perature difference AT (Vt = — 5'qdAT). It can be 
related to the average energy {E) of the carriers by 
5qd = -{E)/eT. 

Transport through a QD depends crucially on the cou- 
pling strength of the dot states to the leads, i.e. the wave 
function overlap of the localized and free states. While 
asymmetries in this coupling have been observed in elec- 
trical transport measurements (see for example [6, 7]), 
their influence on the thermoelectric properties still is 
outstanding. In this paper we compare thermovoltage 
and conductance measurements on a gate-defined QD in 
the Coulomb blockade (CB) regime. In high magnetic 
fields, applied in the plane of the two dimensional elec- 
tron gas (2DEG), a suppression of carrier transport is 
observed for certain single electron tunneling (SET) reso- 
nances. The corresponding thermoelectric signal exhibits 
a strong asymmetry. An analysis of this asymmetry re- 
veals that the QD acts as a thermoelectric rectifier. 

Figure 1 (a) shows a scanning electron microscope im- 
age of the QD structure. This structure is fabricated 
by electron-beam lithography on a GaAs/ (Al,Ga)As het- 
erostructure with a 2DEG carrier density of = 2.3 x 
2q15 ^ nominal mobility of /i = 100 m^/(Vs). 

The QD is embedded at the border of a 2 jiia wide and 20 
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FIG. 1: (a) SEM image of the QD structure; enlarged QD 
region on the right. Schottky-gates are labeled T, R, L, and 
P. (b) Logarithmic grey scale plot (color online) of the con- 
ductance C as a function of B\\ and Vp. (c) Corresponding 
grey scale plot of the thermovoltage (linear scale). Traces at 
0, 7 and 13 T are given separately. 



jiia long electron heating channel [8]. The lithographic 
diameter of the QD is approximately 250 nm and the 
number of electrons can be tuned between 20 and 30 
as a function of applied plunger gate- voltage, Vp. Low- 
frequency {f = 13 Hz) lock-in techniques are used for 
electrical and thermoelectrical measurements. A current 
heating technique is used to create a temperature differ- 
ence of AT ^ 30 mK across the QD. The temperature dif- 
ference and the electron gas temperature (Tg = 80 mK) 
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are determined independently by analyzing SET reso- 
nances as well as the temperature and heating current 
dependence of the universal conductance fluctuations of 
the heating channel. 

Figure 1 (b) presents a grey (color online) scale plot of 
the zero bias conductance in the regime of two SET reso- 
nances ditVp = —2.33 V (resonance A) and Vp = —2.17 V 
(resonance B) as a function of the plunger gate volt- 
age (Vp) and magnetic field {B\\) applied parallel to the 
2DEG plane. Resonance A is typical for standard SET 
CB behavior, whereas the new physics this paper focuses 
on is inferred from the behavior of resonance B. The res- 
onances occur where the chemical potentials of source 
and drain contact leads align with that of the QD, i.e., 
Ms = = Mn+i;n = ^n+1,0 " £^n,o, whcrc E^^i^q and 
En^o denote the ground state (0) energy of n + 1 and n 
electrons on the QD, respectively. Both SET peaks shift 
to less negative gate voltages with increasing magnetic 
fields. This diamagnetic shift is due to a different mag- 
netic field dependence of the energy states in the QD and 
in the leads [6, 9]. More striking is the observation that 
the amplitude of resonance B decreases strongly with in- 
creasing magnetic field and almost vanishes for B\\ > 8 T 
while the amplitude of resonance A remains almost con- 
stant. For clarity, single traces of G{Vp) at = 0, 7, 
and 13 T are shown in the right panel of Fig. 1 (b). 

The corresponding thermovoltage is shown in Fig- 
ure 1 (c). Dark and bright regions correspond to large 
positive and large negative thermovoltage signals, respec- 
tively. SET resonances occur at the sharp transition 
from positive to negative thermovoltages. The observed 
diamagnetic shift is clearly visible in the thermovoltage 
measurement. While the thermovoltage signal of reso- 
nance A shows a uniform behavior for the whole mag- 
netic field range, the thermovoltage of resonance B ex- 
hibits some significant changes for B > 5 T. For this res- 
onance, the characteristic positive- to-negative Vt tran- 
sition turns into a double peak structure with negative 
amplitudes for 5 T < 5 < 8 T. For 5 > 8 T, only a sin- 
gle negative thermovoltage signal remains. In contrast to 
what we observe for this peak in the conductance mea- 
surements, the amplitude of the negative thermovoltage 
signal remains approximately unchanged [right panel of 
Fig. 1 (c)]. In general, such an asymmetric thermovolt- 
age indicates that the electron-hole symmetry of a SET 
process in a QD is broken [10, 11]. 

The broken symmetry is even more obvious in the elec- 
tron addition spectrum, G{Vbs^Vp)^ which is displayed 
in Fig. 2 for 5 = T and 5 = 13 T. In the diamond 
shaped regions, the dot has a constant electron number 
(. . . A/", A/' + l, A/' + 2, . . .). The dark lines (non-zero differ- 
ential conductance) delimiting the diamonds are due to 
the SET resonances, whose positions depend on source- 
drain and plunger-gate voltage. Resonance lines outside 
the CB diamonds originate from transport through ex- 
cited states (/in+i,i;n,j)7 whcrc (i,j)eN+ [12]. The ap- 
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FIG. 2: (color online) Grey scale plot of the differential con- 
ductance as a function of Vbs and Vp at (a) B\\ = T and (b) 
B\\ — 13 T. The arrow marks the transition from a suppressed 
to a high differential conductance. 



pearance of areas with negative differential conductance 
(white areas) is typically a signature of the occurrence 
of blocking mechanisms of various kinds [13, 14]. Espe- 
cially for 5 = 0, the SET resonance /i(7v+2;Ar+i) exhibits 
a negative differential conductance for Fds < V and 
Vp < —2.21 V [onset indicated by an arrow Fig. 2 (a)]. 
The resonance line with a positive slope involving the 
ground state (/iAr+2,0;Ar+i,o), is absent. This blocking of 
the electrical transport for negative Vds indicates that 
the drain contact is aligned with a QD state which cou- 
ples asymmetrically with the contacts. At B\\ = 13 T 
[Fig. 2 (b)], an asymmetric transport gap opens for the 
zero bias + 2 ^ A" + 1 transition, i.e., a new asym- 
metrically coupled state becomes the QD ground state 
A^(Ar+2 Ar+i)' which strongly suppresses the charge carrier 
transport, especially for negative bias voltage. Since in 
our experiments, source and drain contacts are adjusted 
to equal transmission, the observed transport asymme- 
try has to be related to the intrinsic symmetry of the 
QD states. 

The experimental observation can be explained [6, 7] 
by considering two nearly degenerate QD states. The 
5 = ground state, /i(Ar+2;Ar+i), couples symmetrically 
to source and drain contacts while, the second state, 
A^(Ar+2 Ar+i)' couples more strongly to the drain than to 
the source. This situation is schematically depicted in 
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FIG. 3: (color online) SET peak spacings as a function of B\\ 
for constant electron numbers. The curves are distributed 
equally spaced (A = 0.2meV) on a scale between and 
0.6 meV with increasing N. The lines are linear fits with 
slopes of indicated multiples of /is- 



Fig. 4 c). For negative Vds the drain contact will popu- 
late both states and thus the low lying symmetric state 
will become blocked while for positive Vds^ electrons that 
originate from the source will only populate the symmet- 
ric ground state and transport will be possible. 

In a magnetic field parallel to the plane of the 2DEG, 
the Zeeman effect introduces an energy shift to the QD 
states according to Ez = mjgjiBB^ where mj is the 
magnetic quantum number of the QD state with total 
angular momentum J = L -\- g the associated g-factor, 
and jiiB the Bohr magneton. In an in-plane magnetic 
field the orbital component of the electronic wave func- 
tion remains almost unchanged, which ensures that the 
coupling to the leads of the various states is not affected. 
However, the field may very well affect the energetic or- 
dering of the orbital states. From the design of our QD 
one infers that a symmetric coupling to the leads is more 
likely for a symmetric QD state with orbital angular mo- 
mentum L = 0, and that an asymmetric coupling must 
involve states with L 7^ 0. In combination with the above 
observations this implies that, in order to explain the ex- 
periments, a magnetic field must change the energetic 
ordering of these states. 

Information about the energetic ordering of the QD 
states can be obtained by analyzing the magnetic field de- 
pendent energy separation of the SET conductance peaks 
[9, 15]. In Figure 3 the energy spacings between neighbor- 
ing SET resonances is plotted as a function of . The N 
and + 1 states refer to the states with the same index 
in Figs. 1 and 2. For clarity, the curves are distributed 
equally spaced (A = 0.2 meV) on a scale between and 
0.6 meV with increasing N. The curves exhibit a piece- 



wise linear behavior, with slopes that can be expressed 
in integer multiples of gcaAs = (—0-44), i.e., the g-factor 
of the QD states. 

The slope of the SET peak spacings is equivalent to 
the difference in Zeeman energy, AEz = Amjg/jjBB, 
for the two involved QD states. For example, a slope of 
+0.44/i^ is obtained for Amj = mj(^+i) — mj(^) = 1. 
A change of the slope indicates a rearrangement of the 
spin configuration and a state with a different mj*^n-\-i) 
becomes the new ground state. Such rearrangements are 
observed at 5 = 5 and 9.5 T. Of particular interest is 
the increase of the slope for the + 1 state to 0.88 /i^ 
at 5 = 5 T which indicates that indeed a QD state with 
L 7^ is involved in this transitions. 

In order to verify that the simple picture of an asym- 
metrically coupled state can explain the observed ther- 
movoltage signal, we have calculate the electric, ther- 
moelectric and thermal transport coefficients within the 
Laudauer formalism of resonant tunneling [16]. The total 
generalized current through the quantum dot is given by 



tot 



= fdE( 

J —CX) \ 



^]{fL{E,T) 



ME,T)]t{E) (1) 



where /l (fn) is the Fermi distribution function of the 
left (right) contact, t{E) is the energy dependent trans- 
mission coefficient of the QD, and A is either the charge 
(— e) or the energy {E — fi) which is transported through 
the QD, as applicable to for charge or heat currents, re- 
spectively. The charge (/) and heat (Q) currents are 
related to the applied electrochemical potential and tem- 
perature difference via the transport coefficients Lij : 
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where Ln = G and the thermopower is defined by the 
quotient S = —{L12/L11). 

The energy dependent transmission coefficient can be 
modeled as follows: 



t{E,T) = A 



(r/2)2 



(r/2)2 + (^)2 



f{E-6Ez^T), (3) 



The first term of Eq. 3 represents a single QD resonance, 
while the second term accounts for the (thermal) occu- 
pation of the blocking state which depends on the energy 
separation 5Ez between the transmitting and the block- 
ing state. 

Figure 4 (a) shows the calculated conductance for 
various energy separations. The values of A = 0.79 
and F = 0.2 meV result from a fit of resonance B for 
= T. A constant background Gcot = 0.001 e'^/h has 
been added in order to account for co-tunneling contribu- 
tions via additional excited states. The amplitude of Gcot 
corresponds to the experimentally observed minimum 
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FIG. 4: (color online) Calculated conductance G (a) and ther- 
mopower S (b) for Te = 80 mK, AT = 30 mK and vari- 
ous energy separations between the symmetric and asymmet- 
ric (blocking) states: (1) without blocking state, (2) 3Ez = 
—250 /ieV, (3) /xeV, and (4) 400 /xeV. The crosses represent 
the measurement B\\ — 13 T. (c) presents a schematic en- 
ergy diagram of the QD close to the + 1 ^ + 2 SET 
resonance at By = T and 13 T. The blocking state is sym- 
bolized by the thick broken line, (d) displays the difference 
of the thermal conductances A^t for a reversal of the temper- 
ature difference. 

conductance between two SET conductance peaks. Fur- 
thermore, the experimentally determined electron tem- 
perature, Tg = 80 mK, and temperature difference of 
AT = 30 mK have been used for all calculations. 

Starting from a situation where the blocking state 
is energetically above the symmetric state {5Ez = 
—250 jieV) the SET conductance peak becomes more and 
more asymmetric with increasing 5E (c.f. Fig. 4 (a)). For 
5Ez = a sharp suppression of the conductance at the 
center of the SET peak is observed. The calculated con- 
ductance for 5Ez = 400 /ieV [curve 4 in Fig. 4 (a)] closely 
resembles the measured G{Vp) at 5|| = 13 T. This value 
of SEz corresponds to a situation where the blocking 
state is energetically below the symmetric state [Fig. 4 
(c)]. The corresponding thermopower [Fig. 4(b)] changes 
from a resonance-like line shape to a single peak struc- 
ture, where good agreement with the experiment (dashed 
line) is again achieved for SEz = 400 jaeV. 

The change of the line shape for the thermopower is 
direct evidence that the QD acts as a thermal recti- 
fier. Fig. 4(d) shows the difference between the ther- 
mal conductance, = /^TiTs — ^TsTi, with = L22 — 
{L12L21/ Lii) and the Lij obtained from the fits in Figs. 
4 (a) and (b), for a reversal of the temperature difference 
across the QD and SEz = 400 fieV. At Vp = -2.11 V 
the efficiency of rectification, {Ak/k:) reaches 10.5% at 



the given temperatures. This number, while modest, is 
significantly higher than the rectification achieved in Ref. 
[4], and can be even further increased by an optimized 
design of the QD layout. 

In summary, we have demonstrated that in a QD asym- 
metric coupling to its leads not only modifies the charge 
transport, but also strongly modifies its thermal proper- 
ties. If the QD design is such that transport is favorable 
through states with non-zero orbital momentum, high ef- 
fective thermal rectifying properties can be obtained. In 
contrast to the previously reported nanotube based ther- 
mal rectifiers [4], QDs directly control the heat transfer 
in the electronic system, without the need for additional 
coupling to the phonon system. This opens the perspec- 
tive of more sophisticated electronic devices with high 
rectifying performances. 
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